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The displacement- and frequency-noise-free interferometer (DFI) is a multiple laser interferometer 
array for gravitational wave detection free from both the displacement noise of optics and laser 
frequency noise. So far, partial experimental demonstrations of DFI have been done in 2-D table top 
experiments. In this paper, we report the complete demonstration of a 3-D DFI. The DFI consists 
of four Mach-Zehnder interferometers with four mirrors and two beamsplitters. The displacement 
noises both of mirrors and beamsplitters were suppressed by up to 40 dB. The non- vanishing DFI 
response to a gravitational wave was successfully confirmed using multiple electro-optic modulators 
and computing methods. 



Gravitational waves are perturbations of the space- 
time curvature propagating across the universe at the 
speed of light. They propagate in space-time, extending 
a proper distance between test masses in one transverse 
direction, and shortening in the other, orthogonal direc- 
tion. A laser-interferometric gravitational wave detector 
measures length differences between the optical compo- 
nents (test masses) caused by the gravitational wave ef- 
fect. Various noise sources that limit the sensitivity of the 
detector must be removed to achieve the best sensitivity 
level, as the interferometers must detect the very tiny 
signals of length change. The sensitivities of the ground- 
based detectors are affected by seismic noise, thermal 
noise, radiation pressure noise, laser source noise (laser 
frequency and intensity noise), and shot noise. The seis- 
mic, thermal and radiation pressure noises can be catego- 
rized as displacement noises since the noise sources shake 
up the optics directly. These displacement noises limit 
the sensitivities especially in the low-frequency regime, 
which is the important region for the ground-based grav- 
itational wave detectors, which target at several ten to 
several thousand Hz. Recently, the idea of a gravitational 
wave detector free from both the displacement noises of 
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the optics and the laser frequency noises at all frequen- 
cies has been proposed in Refs. p], Q. This technique 
exploits the fact that the gravitational wave effect on the 
laser light takes a form different from that of optics dis- 
placements. A signal combination which does not sense 
displacement- or frequency-noises but senses the gravi- 
tational wave contribution can be constructed when the 
array of N test masses consists of multiple interferome- 
ters under the condition that N > d + 2 (d is the number 
of spatial dimensions of the array). The next paper in the 
series of DFI papers, Ref. 0, suggested a 3-D DFI con- 
figuration consisting of four Mach-Zehnder interferome- 
ters (MZIs). This configuration is the most ideal of all 
DFI configurations so far proposed. A proof-of-principle 
experiment with a 2-D layout of a sub-part of this con- 
figuration was done in Ref. [4|]. There, the cancellation 
of displacement noise from a folding mirror was demon- 
strated. However, the 2-D layout was half that of the 3-D 
DFI, and it cannot cancel the displacements of all the op- 
tics. In this paper, as a complete DFI system, we report 
the cancellation of the displacement both of the folding 
mirrors and beamsplitters while the gravitational-wave 
response is retained using the full 3-D DFI. 

The 3-D DFI configuration is depicted in Fig. [TJ Four 
MZIs with equal arm lengths compose the DFI. Let us de- 
note by IFOi, IF0 2 , IF0 3 , and IF0 4 the four single MZIs 
drawn in orange (ACiB- ADxB), pink (AC 2 B- AD 2 B), 
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FIG. 1: (color online). Optical topology of 3-D DFI. The 
pairs of IFOi and IFO3 and of IFO2 and IFO4 are coun- 
terpropagating to each other, constructing two bidirectional 
Mach-Zehnder interferometers (MZIs). Right and left bidi- 
rectional MZIs share the two beamsplitters. The length from 
A or B to each folding mirror is L. 



light blue (BdA - BD X A), and green (BC 2 A - BD 2 A) 
in Fig. [H Each length between a beamsplitter and a mir- 
ror (AC 1 ,AC 2 ,AD 1 ,AD 2 ,C 1 B,C 2 B,D 1 B, and D 2 B) is 
assumed as L so that each interferometer is insensitive to 
laser frequency noise thanks to the equal arms. Two of 
the four interferometers are combined to construct one 
bidirectional MZI in which two beams are counterpropa- 
gating on the same path. As is discussed in Ref. [3|], the 
DFI signal combination can be written as 

0DFI = (01 - 03) - (02 - 04) (1) 

where 0j(j = 1, 2, 3, 4) are the phase signals of IFOj. In 
Eq. (pQ), the combined signal of (0i — 03) is the output 
of the bidirectional MZI that is free from the displace- 
ments of C\ and D\. Because the folding mirrors are 
at the midpoint of the arms, and the laser beams travel 
in both directions in the MZI, the phase variation due 
to the displacement of C\ and D\ arrives at each output 
port simultaneously. Therefore, in the frequency domain, 
IFOi and IFO3 have identical transfer functions from the 
displacement of C\ and D\ to the signal ports, expressed 
as 

H? liDl (Q) = H^' Dl (Q) = e~ iLn/c , (2) 

where c is the speed of light in vacuum, uo is the laser fre- 
quency and fl the Fourier frequency of the displacement 
motion. The phase modulations due to the motions of the 
optics experience the same phase delay shown in Eq. ([2]) 
because the displacement sources are at the midpoint of 
the path. These identical responses to C\ and D\ can 
thus be removed by subtraction. Similarly, (02 — 04) is 
the signal of the other pair of bidirectional MZIs, and 
free from the displacements of C 2 and D 2 . Any signals 
that phase modulate the light and do not arise at the 
exact center of the path will not arrive simultaneously at 
the outputs, and thus will not be canceled by the bidi- 



rectional MZIs. The beamsplitter displacements are re- 
moved by the combination of (0i — 4> 2 ) in Eq. (pQ). Since 
the laser beams of IFOi and IFO2 reach the two beam- 
splitters simultaneously, their responses to the beamsplit- 
ter motion are the same. Therefore, in the frequency do- 
main, IFOi and IFO2 have identical transfer functions 
from the displacement of A to the signal ports expressed 
as 

Hf(Q) = H£(Q) = e~ 2iLQ / c . (3) 

The transfer functions of the displacement motion of 
beamsplitter B are i/f (ft) = H^(Q) = 1 since the mo- 
tions of B are just before the detection ports. Therefore 
the beamsplitter displacements can be eliminated by the 
combination of IFOi and IFO2. Similarly, the combined 
signal of (—03 + 04) in Eq. (pQ) is free from the beamsplit- 
ter displacements. In practice, in addition to the trans- 
fer functions, extra phase delays should be considered 
that arise from the paths from the second beamsplitter 
to the signal analyzer via a photo detector and signal 
cables. These additional phase delays of each interfer- 
ometer should be made equal or calibrated later, so as 
not to lose the balanced responses. 

According to Ref. [3], the DFI signal responds to a 
gravitational wave with the transfer function, 

ffcw(tt) = " 77 j (! + 1/V2)e-?Z (1 - e 

-(1 - l/V2)eVTc (1 - e L ± J — L ) I 

where ft is the Fourier frequency of the gravitational wave 
and we assume the gravitational wave of 77 — £ polariza- 
tion is coming along the z direction in Fig. [U The 3-D 
DFI does not respond to DC gravitational waves because 
the response is degraded by the subtraction when the 
DFI signal is built by the four MZI signals. The DFI 
response is proportional to (ftL/c) 2 below the peak fre- 
quency. The peak frequency depends on the arm length 
L, e.g. when L = 1 m, the peak is at about 150 MHz. 
There have been several attempts to lower the effective 
frequency of DFI: Adding Fabry-Perot (FP) cavities to 
the 3-D DFI to increase the effective path lengths and 
lower the peak frequency; another method of displace- 
ment cancellation using one FP cavity @, 0- At this 
moment, our f 2 response at low frequencies is the best 
sensitivity that DFI configurations have achieved. 

Figure. [21 shows the practical setup of our experiment. 
The laser source is a commercial solid-state Nd:YAG laser 
at 1064 nm. The light source is split into four by beam- 
splitters so as to provide the four incident beams for the 
four MZIs. As depicted in Fig. El IF0 2 and IF0 4 (de- 
noted as higher interferometers) are lying over IFOi and 
IFO3 (denoted as lower interferometers). The length L 
was chosen to be about 40 cm yielding peak frequency 
of about 370 MHz. There are two Faraday isolators be- 
fore the lasers are sent to IFOi and IFO2. They are 
used for extracting the signal of IFO3 and IFO4 (ports 




FIG. 2: (color online). Schematic of the experiment. The 
octahedron is pushed over sideways to be set up on an optical 
table. IFOi and IFO3 comprise one bidirectional MZI, and so 
do IFO2 and IFO4. EOM1 simulates mirror C\ displacement 
and EOM2 and EOM3 simulate beamsplitter B displacement. 




FIG. 3: (color online). The overhead view of the experimental 
setup. HWP: half wave plate, QWP: quarter wave plate, FI: 
Faraday islator, BS: beamsplitter, DCPD: DC photo detector. 
There are four interferometers: IFOi (A — Ci,Di,—B), IFO2 
(A - C 2 ,D 2 ,-B), IFO3 (B - Ci,L>i,-C) and IF0 4 (B - 
C2, -D2, — A). The black-dashed line indicates the optical lever 
method for the alignment of 3-D configuration. 

U3 and U4 [8[ exiting from the IFOi and IFO2 inputs, 
respectively. The signals at ports trl and tr2 are directly 
detected by the photo detectors. The DC photo detectors 
(DCPD) are placed at ports trS and trA for extracting 
the control signals. The extracted DC signals are sent 
to the servo filters with static offsets, then sent to the 
piezoelectric transducers (PZTs) attached to mirrors D\ 
and C2. D\ and C2 are controlled by signals giving the 
mid- fringe locking for the IFO3 and IFO4, respectively. 
Once IFO3 is locked, IFOi is locked automatically. Like- 
wise, locking IFO2 automatically locks IFO4 The control 
bandwidth of the two interferometers is about 1 kHz. 

To align the four-interferometer system, we utilized an 
optical lever as shown in Fig. 03 using the fact that the 
lower and higher sets of bidirectional MZIs share their 
beamsplitters [7[ . The positions and the angles of C\ and 
of Di, and positions and angles of C2 and of D2, inde- 
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FIG. 4: (color online). The solid green line is the transfer 
function from EOM3 to tt4 as a reference, and the solid black 
line is the transfer function from EOM3 to the DFI signal. 
From one to two orders of magnitude of the cancellation were 
attained in the DFI signal. 



pendently determine the position and angle of B for the 
alignment of the higher and lower interferometer, respec- 
tively. To make them compatible, the optical lever was 
utilized to monitor the angle of beamsplitter B. First, 
for the lower interferometer, the angle of B determined 
by C\ and D\, is adjusted and monitored by the optical 
lever signal on the screen For the higher interferometer, 
we then adjust the position and angle of B determined 
by C2 and D2, is adjusted for the lower interferometer 
and displayed on the screen. And we adjust the posi- 
tions and angles of C2, D\ and A so that the two moni- 
tored marks on the screen come to exactly the same point 
which means the ideal angle of B. This method achieved 
an almost perfect contrast. 

EOMs were utilized to create phase shifts to the laser 
light as simulated displacements of the optics because the 
interesting frequency, i.e., the response peak of the grav- 
itational wave, is expected to be in the hundred-MHz re- 
gion. To confirm the noise cancellation features of both 
mirror and beamsplitter, the three EOMs were driven; 
for simulating the displacement of mirror Ci, EOM1 is 
placed at the center of the path BC\A\ for simulating 
the displacements of beamsplitter B, EOM2 and EOM3 
are placed close to B as beamsplitter B affects both the 
higher and lower interferometer. The EOMs were driven 
by swept-sine noise sources provided by the RF network 
analyzer (Agilent 4395 A). The signals from the four MZIs 
are detected by four fast photo detectors (New Focus 
1611) at ports £rl,tr2,££3 and ttA. The optical gain im- 
balances between the four output signals are compen- 
sated by neutral density (ND) filters placed in front of 
each detector. The DFI signal is obtained by electrical 
summation of the four signals using electric power com- 
biners. We adjusted the control polarity so that trl and 
tr2 have opposite sign so as to cancel by their summation. 

The DFI response to a gravitational wave coming along 



4 




io 6 io 7 

Frequency [Hz] 



FIG. 5: (color online). The orange dots are the simulated 
DFI response to the gravitational wave obtained by being 
measured and computed, the gray solid line is the DFI re- 
sponse predicted by Eq. Jl}. The simulated response agreed 
with the theory. 

z direction was demonstrated by multiple measurements 
and simulation. Ideally, the continuous phase variation 
should be applied on the entirety of the laser paths to 
simulate the gravitational- wave effect. However, only the 
part covered by the crystal causes a phase shift by an 
EOM. Therefore, the following procedure was adopted. 
We place an EOM in the laser path and measure the 
transfer function. Then we place the EOM at the next 
point and measure the transfer function, and so on. The 
measured data are summed up with proper phase delays 
that depend on the gravitational wave direction and the 
EOM position. For a gravitational wave coming along the 
AB axis, Eq. (|4]) can be approximated by the summation 
of the discrete phase shifts expressed as 

H GW (Q) ~ H sum (Q) 

T , 2D 2D 

(jO±j I iQL ^ — > jflL n_ _i£ L L L n_ iflL ^ — > iflL n_ i^L n 

— 1 g \/2c \ g c D g \/2c D — g \/2c \ g c D g \/2c D 



under the condition QL/Dc < 1. D is the resolution 
of the positions on the path of L where the phase shifts 
occur. Note that the straightforward way to produce 



continuous phase shift is to insert EOMs along all laser 
path. However, that many EOMs covering laser paths 
will cause a serious reduction in interferometer contrast. 

In the experiment, the phase shift part of e lQLn / Dc 
in Eq. (O was created practically by the EOM with 
the discrete resolution D = 10. For the first ten data, 
n = ... 9, an EOM was set in a respective position at 
nL/D from beamsplitter £?, and the respective trans- 
fer functions were measured. Since we cannot put an 
EOM on C\A, we compute the transfer function data for 
n = 10 ... 20 applying proper phase delays to the data of 
n = 9. These data were summed with predicted phase 
delays due to the incident angle of the gravitational wave 
in accordance with Eq. ([5]), to construct the gravitational 
wave response. 

The results for displacement-noise suppression are 
shown in Fig. [H The solid black line is the transfer func- 
tion from EOM1 to the DFI (combined) signal. The solid 
green line is the transfer function from EOM1 to a single 
MZI output (port UA) as a reference. Comparing with 
these transfer functions, we can confirm that the DFI sig- 
nal has less sensitivity to the displacement by one to two 
orders of magnitude. The transfer functions in Fig. [H 
include the response of the EOMs and the photo detec- 
tors, and also phase delays due to the optical and signal 
paths after the photo detectors to the analyzer. 

In Fig. [5j the orange dots are the simulated DFI re- 
sponse to a gravitational wave obtained by the measure- 
ment and computing, and the gray solid line is the DFI 
response predicted by Eq. The response function of 
photo detectors, and of the mechanical resonance due to 
the EOM crystal, and the phase delays due to the op- 
tical paths and signal cables from beamsplitters to the 
analyzer are compensated properly in Fig. [5j The 3- 
D DFI response to a gravitational wave is well demon- 
strated around the peak frequency and its f 2 feature can 
be confirmed just below that frequency. 

To summarize the results: we have successfully con- 
structed and operated the 3-D DFI. The noise suppres- 
sion feature of about 1 to 2 orders of magnitude and the 
retained DFI response were well simulated using EOMs 
over a wide frequency range. 
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